home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / SPLIN2.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  669b  |  23 lines

  1. PROCEDURE splin2(x1a,x2a: glnarray; ya,y2a: glmbyn;
  2.        m,n: integer; x1,x2: real; VAR y: real);
  3. (* Programs using routine SPLIN2 must define the types
  4. TYPE
  5.    glnarray = ARRAY [1..nn] OF real;
  6.    glmbyn = ARRAY [1..m,1..n] OF real;
  7. in the main routine. The dimension nn of glnarray must be set
  8. to the larger of n and m *)
  9. VAR
  10.    k,j: integer;
  11.    ytmp,y2tmp,yytmp: glnarray;
  12. BEGIN
  13.    FOR j := 1 TO m DO BEGIN
  14.       FOR k := 1 TO n DO BEGIN
  15.          ytmp[k] := ya[j,k];
  16.          y2tmp[k] := y2a[j,k]
  17.       END;
  18.       splint(x2a,ytmp,y2tmp,n,x2,yytmp[j])
  19.    END;
  20.    spline(x1a,yytmp,m,1.0e30,1.0e30,y2tmp);
  21.    splint(x1a,yytmp,y2tmp,m,x1,y)
  22. END;
  23.